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Abstract 

^ Loop percolation, also known as the dense O(l) loop model, is 

a variant of critical bond percolation in the square lattice I? whose 
graph structure consists of a disjoint union of cycles. We study its 
connectivity pattern, which is a random noncrossing matching as- 
sociated with a loop percolation configuration. These connectivity 
patterns exhibit a striking rationality property whereby probabilities 
■^j- of naturally-occurring events are dyadic rational numbers or rational 

functions of a size parameter n, but the reasons for this are not com- 
pletely understood. We prove the rationality phenomenon in a few 
cases and prove an explicit formula expressing the probabilities in the 
^ "cylindrical geometry" as coefficients in certain multivariate polyno- 

t-h mials. This reduces the rationality problem in the general case to 

that of proving a family of conjectural constant term identities gener- 
alizing an identity due to Di Francesco and Zinn-Justin. Our results 
make use of, and extend, algebraic techniques related to the quantum 
Knizhnik-Zamolodchikov equation. 
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1 Introduction 



1.1 Critical bond percolation and loop percolation 

Critical bond percolation on Z 2 is the most natural model of a random sub- 
graph of the square lattice: each edge of the lattice is included with prob- 
ability 1/2, independently of all other edges. Despite the simplicity of its 
definition, it is difficult to analyze rigorously. Though the spectacular recent 
results of Smirnov [26J and Lawler- Schramm- Werner [17] related to critical 
exponents, conformal invariance and the Schramm-Loewner Evolution (SLE) 
have resolved many of the outstanding open problems for the related model 
of critical site percolation on the triangular lattice, the same problems are 
still unsolved in the case of bond percolation on Z 2 . 

In this paper we will study another percolation model, which is related to 
critical bond percolation on Z 2 and can be thought of as a variant of it. We 
will refer to the model as loop percolation on Z 2 ; an equivalent model has 
been studied in the statistical physics literature under several different names 
such as the dense 0(1) loop model [TH], completely packed loops [33] 
and the XXZ spin chain [3j E] . 

Our motivation in studying this model was twofold: first, to point out, 
and attempt to systematically exploit, the inherent interest and approach- 
ability of the model from the point of view of probability theory — qualities 
which have not been emphasized in previous studies. Second, to use and en- 
hance some of the deep algebraic tools (referred to under the broad heading 
of the quantum Knizhnik-Zamolodchikov equation) that were devel- 
oped in recent years to compute explicit formulas for several quantities of 
interest in the model. In this paper we focus on algebraic techniques. The 
follow-up paper [21] will contain additional results of a more probabilistic 
flavor. 

Let us start with the definition of the model. We define loop percolation 
in the entire plane, but later we will also consider the model on a half-plane 
and on a semi-infinite cylinder. 

Definition 1.1 (Loop percolation). Consider the square lattice Z 2 equipped 
with a checkerboard coloring of the squares of the dual lattice, such that the 
square whose bottom-left corner is (0, 0) is white. Loop percolation is the 
random subgraph LP = (Z 2 ,V^(LP)) of Z 2 obtained by tossing a fair coin 
for each white square □ in the dual lattice, independently of all other coins, 
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and including in V(LP) either the bottom and top edges of □ or the left and 
right edges of □ according to the result of the coin toss (Fig. [if. 



• 9 * • 

Figure 1: The two possible edge configurations around a white lattice square. 

Fig. [2] shows portions of a critical bond percolation configuration and a 
loop percolation configuration. In fact, the two models are equivalent, in 
the following sense. The black squares form the vertices of a graph where 
two black squares are adjacent if their centers differ by one of the vectors 
(0, +2), (+2, 0). Identifying each black square with its center, this graph 
decomposes into disjoint components of "blue" and "red" sites according to 
parity, each component being isomorphic to Z 2 . Define a subgraph G of the 
blue component whose edges are precisely those that do not cross an edge 
of the loop percolation graph. Then it is easy to see that G is a critical 
bond percolation graph; see Fig. [2j Conversely, it is clear that starting from 
the critical bond percolation graph one can reconstruct the associated loop 
percolation graph. 

As an immediate consequence of the above discussion, we see that the 
loop percolation graph almost surely decomposes into a disjoint union of 
cycles, or "loops." Indeed, since each site of Z 2 is incident to precisely two 
white squares, it has degree 2, so the connected components of the graph are 
either loops or infinite paths; however, the presence of an infinite path would 
imply the existence of an infinite connected component in the "blue" critical 
bond percolation graph containing those blue sites adjacent to an edge in the 
infinite path, in contradiction to the well-known fact [T3l Lemma 11.12] that 
critical bond percolation on Z 2 almost surely has no such infinite component. 

Given the equivalence between loop percolation and critical bond per- 
colation described above, one may ask why it is necessary to study loop 
percolation separately from critical bond percolation. We believe that there 
is much to be gained in doing so. In particular, the existing research on loop 
percolation and connections with other natural statistical physics models 
such as the XXZ spin chain and Fully Packed Loops indicate that studying 
this type of percolation as a model in its own right suggests a variety of 
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Figure 2: (a) Critical bond percolation, (b) Loop percolation. 
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Figure 3: The critical bond percolation graph associated with a loop perco- 
lation graph. 



mathematically interesting questions that one might not otherwise be led to 
consider by thinking directly of critical bond percolation. 

Furthermore, the techniques developed for studying loop percolation have 
shown that it belongs to the (loosely defined) class of so-called "exactly solv- 
able" models, for which certain remarkable algebraic properties can be used 
to get precise formulas for various quantities associated with the model. The 
new results presented in this paper will provide another strong illustration of 
this phenomenon and give further credence to the notion that loop percola- 
tion is quite worthy of independent study. Ultimately, we hope that this line 
of investigation may lead to new insights that could be used to attack some 
of the important open problems concerning critical bond percolation. 

1.2 Noncrossing matchings 

Our discussion will focus on certain combinatorial objects known as non- 
crossing matchings that are associated with loop percolation configura- 
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tions. Let us recall the relevant definitions. For two integers a < b, let [a, b] 
denote the discrete interval {a, a + l,...,b}. Recall that a noncrossing 
matching of order n is a perfect matching of the numbers 1, . . . , 2n (which 
we encode formally as a function vr : [1, 2n] — > [1, 2n] such that ttott = id and 
ir(k) k for all k e [1, 2n]; ir(k) represents the number matched to k), which 
has the additional property that there do not exist numbers a < b < c < d 
in [1, 2n] such that it (a) = c, n(b) = d. If ir(j) = k we say that j and k 
are matched under tt and denote j * — > k. Denote by NC n the set of 
noncrossing matchings of order n. It is well-known that the number |NC n | 
of elements in NC n is Cat(n) = ^y( 2 ^), the nth Catalan number. 

Similarly, an infinite noncrossing matching is a one-to-one and onto 
function tt : Z — > Z that satisfies tt o tt = id and ir(k) k for all k e Z, such 
that there do not exist integers a < b < c < d for which tt matches a to c 
and b to d. Denote by NC^ the set of noncrossing matchings on Z. 

Both finite and infinite noncrossing matchings can be represented graph- 
ically in a diagram in which noncrossing edges, or arcs, are drawn between 
the elements of any matched pair j *—> k; see Fig. (ij As illustrated in the 
figure, in the case of a finite noncrossing matchings there are two equivalent 
representations, as a matching of 2n points arranged on a line or on a circle. 

1.3 The connectivity pattern of loop percolation: the 
infinite case 

We now associate a random noncrossing matching called the connectivity 
pattern with loop percolation configurations. There will be two variants 
of the problem according to the type of region being considered, resulting 
in infinite or finite noncrossing matchings. We start the discussion with the 
infinite case. To the best of our knowledge, this variant of the problem has not 
been previously considered in the literature, but nonetheless the rationality 
phenomenon we will discuss appears most striking in this setting. 

Definition 1.2 (Half-planar loop percolation). Denote by Z^, E the half- 
lattice Z^ E = {(m, n) : m + n ^ 0}, and consider a loop percolation graph 
LPne defined as before but only on Z^ E instead of on the entire plane. The 
connectivity pattern associated with it is an infinite noncrossing matching, 
which we denote by IT*. It is defined as follows: for any j e Z, to compute 
n*(j), start at the vertex (j, — j). According to our convention regarding the 
coloring of the squares of the dual lattice, (j, — j) is at the bottom-left corner 



6 



9 8 




-5 -4 -3 -2 -1 1 2 3 4 5 6 7 8 9 10 
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Figure 4: (a-b) a finite noncrossing matching shown as a matching of points 
on a line or on a circle; (c) an infinite noncrossing matching. 

of a white square, so (j, —j) is incident to precisely one edge in the configu- 
ration. Follow that edge, and, continuing along the path of loop percolation 
edges leading out from one eventually ends up at a vertex (k,—k) 

on the boundary of the half-lattice (since if we consider the configuration as 
part of a configuration on the entire plane, the path must be a loop as noted 

above). In this case, we say that j and k are matched, and denote j k; 
see Fig. [5] for an illustration. 

It is clear from elementary topological considerations that IT* is a (ran- 
dom) noncrossing matching. It turns out to have some remarkable distri- 
butional properties. In particular, a main motivation for the current work 
was the observation that the probabilities of many events associated with II* 
turn out rather unexpectedly to be rational numbers which can be computed 
explicitly. We refer to this as the rationality phenomenon. We will prove 
it rigorously in a few cases, and conjecture it for a large family of events. 
Although we have not been able to prove the conjecture in this generality, 
we also derive several results that establish strong empirical and theoretical 
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(b) 

Figure 5: (a) A loop percolation configuration (the paths leading out of the 
vertices (n, —n) are highlighted and colored to emphasize the connectivities); 
(b) the associated connectivity pattern. 
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evidence for the correctness of the conjecture and reduce it to a much more 
explicit algebraic conjecture concerning the Taylor coefficients of a certain 
family of multivariate polynomials. 

Given a finite noncrossing matching n e NC n and an infinite noncrossing 
matching 7r 6 NC^, we say that 7r is a submatching of 7r if 7Ti = 7r , 

|[l,2n] 

and in this case denote tt < ir. We refer to an event of the form {ti < II*} 
as a submatching event. We will often represent a submatching event 
schematically by drawing the diagram associated with the submatching ttq] 
for example, the diagram "-^V-^V" corresponds to the submatching event 
{(n*) [M] = (2,1,4,3)}. 

Conjecture 1.3 (The rationality phenomenon for submatching events). For 
any 7r e NC n; the probability P(vr < II*) of the associated submatching 
event is a dyadic rational number that is computable by an explicit algorithm 
(Algorithm C described in Appendix A). 

Table [T] lists some of the simplest submatching events and their conjec- 
tured probabilities. We can prove the first two cases. 

Theorem 1.4. We have 

(YY« n -H- « 
< n.) - ^ (2) 

As we explain in Subsection 1.7 the relation Q will follow as an imme- 

in 



P 




diate corollary of a result due to Fonseca and Zinn- Justin (Theorem 1.11 



Subsection 1.5). The second relation (|2| can also be derived from the same 
theorem but in a slightly less trivial manner; see Theorem 1.16 in Subsec- 
tion 11.61 

We note that the rationality phenomenon is not restricted to submatch- 



ing events; a similar statement to Conjecture L3 also appears to hold for a 
larger family of finite connectivity events, which are events of the form 

k)eA \j for some finite set A a Z x Z. For example, in Theo- 



rem IT7T8 in Subsection 1.7 we derive rigorously the value 135/1024 as the 



probability of a certain finite connectivity event that is not a submatching 
event. However, we mostly focus on submatching events since our main re- 



sults (Theorems 1.13 and 2.24) pertain to them and as a result provide strong 



theoretical evidence for a rationality phenomenon in this setting. 
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Event 


Probability 




Q 
O 


1 2 


8 




97 


12 3 4 


512 




59 


12 3 4 


1024 



Event 



^ ^ f\ 

1 2 3 4 5 6 



1 2 3 4 5 6 



1 2 3 4 5 6 



1 2 3 4 5 6 




Probability 



214093 

2 21 

69693 
2 21 

69693 
2 21 

37893 
2 21 

7737 



2 21 



Table 1: The probabilities of some submatching events in the noncrossing 



matching II*. The first two are established rigorously (Theorem 1.4), others 
are conjectured. 

1.4 Loop percolation on a cylinder: background 

Our study of half-planar loop percolation and its connectivity pattern is based 
on the fact that the model can be approached in a fairly straightforward 
manner as a limit of a model with a different geometry of a "semi-infinite 
cylinder." This corresponds to loop percolation in a semi-infinite diagonal 
strip whose two infinite boundary edges are identified. The precise definition 
is as follows. 

Definition 1.5 (Cylindrical loop percolation). For any n ^ 1, the cylindrical 
loop percolation graph LP n has vertex set 



{(x, t/)eZ 2 : x + y ^ 0, —n + 1 ^ x — y < n}, 



(3) 



where for any k 5= we identify the two vertices (— n + 1 + k,n — 1 + k) and 
(n + k, —n + k). The edges are sampled randomly in the same manner as the 
usual loop percolation. 
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There is a convenient way of representing the cylindrical model as a tiling 
of a strip of the form [0, 2n] x [0, go) with two kinds of square tiles, known as 
plaquettes, which are shown in Fig. [6j To see the correspondence, first ro- 
tate the diagonal strip in the original loop percolation configuration counter- 
clockwise by 45 degrees, then replace each (rotated) white lattice square with 
a plaquette whose connectivities correspond to the percolation edges in the 
white square. Once the plaquette tiling representation is drawn, one can 
wrap it around a cylinder to obtain a three-dimensional picture; see Fig. [7j 
The representation using plaquettes is the one traditionally used in most of 
the existing literature. 




Figure 6: The two types of plaquetfr 



Let H* denote the connectivity pattern associated with the cylindrical 
loop percolation graph LP n , defined analogously to H* by following each of 
the paths originating at the vertices {—j,j),—n + 1 < j < n until it re- 
emerges in a vertex (—k,k). It is easy to see that II* is a random finite 
noncrossing matching of order n; the only justification needed is the following 
trivial lemma that shows that it is well-defined. 

Lemma 1.6. In the graph LP n , almost surely all paths are finite. 

Proof. Consider the tiling representation of the model as in Fig. |7j^b). It is 
easy to see that a row in the tiling in which the plaquettes alternate between 
the two types of plaquettes forces all paths below it to be finite by "bouncing 
back" any path attempting to cross the row. Any row has probability l/2 2n ~ 1 
to have such structure, independently of other rows, so almost surely there 
will be infinitely many such rows. □ 

The cylindrical connectivity pattern IT*™'' has been the subject of extensive 
research in recent years, and appears in a surprising number of ways that 
do not seem immediately related to each other. Before presenting our new 
results on the behavior of let us survey some of the known theory. 
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Figure 7: (a) Cylindrical loop percolation: the north-west and south-east 
boundary edges of the infinite diagonal strip are identified along the dashed 
lines; (b) representation of the same configuration as a tiling of two kinds of 
plaquettes; (c) wrapping a tiling of plaquettes around a cylinder. 
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For a matching 7r e NC n denote fi n = = it). A natural ques- 

tion is how to compute the probability vector fi n = (/i 7r ) 7r eNC„- It is easy to 
see that \i n is the stationary distribution of a Markov chain on NC„ where 
the transitions n — » n' correspond to the operation of extending the cylin- 
der by one additional row of random plaquettes. This operation leaves the 
distribution fi n invariant since it results in a plaquette tiling equal in dis- 
tribution to the original one. Formally, there is a Markov transition matrix 
(^C^WeNCn such that 



J- n 



(1/2) 



A^n^n ^ ^ Mm! 



where fi n is considered as a row vector. (In the statistical physics literature 
Tn^ is usually referred to as a transfer matrix.) 

The reason for the notation Tn^ is that we can generalize the model 
and define a matrix for any < p < 1, corresponding to a tiling of 
independently sampled random plaquettes in which each choice between the 
two types of plaquette is made by tossing a coin with bias p. Remarkably, 
the value of p does not affect the distribution of the connectivity pattern. 

Theorem 1.7. The transition matrices (T^) 0<p< i are a commuting family 
of matrices. Consequently, since they are all stochastic, they all share the 
same row eigenvector \x n associated with the eigenvalue 1. 

Note that using the transition matrix Tn^ 2 \ or for any fixed p, is 
not necessarily the easiest way to compute /u n , since to compute the entries 
of Ti 1/2) one has to count the number of possible rows of n plaquettes (out 
of the 2 2n possibilities) that would cause a given state transition 7r — > n'. It 
turns out that there is a more convenient (and more theoretically tractable) 
system of linear equations for computing fi n , which involves a simpler matrix 
H n that arises as a limiting case of the ri p) as p 
p — > 1. 

To define H n , first define for each k e [1, 2n] a mapping : NC 
given by 



0, or symmetrically as 



e fc (7r)(m) 



n{m) 


if m 


k + l 


if m 


< k 


if m 


n(k + 1) 


if m 


n(k) 


if m 



k + 1, 

7T(k), 



(4) 



1), 
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where k + 1 is interpreted as 1 if k = 2n. In words, efc(7r) is the matching tt' 
obtained from ir by unmatching the pairs k <— > ir(k) and k + 1 <-^-> 7r(A; + 1) 

and replacing them with the matched pairs *—*■ k+1 and 7r(/c) <-^— >• 7r(/c+l). 
In the case when and + 1 are already matched under 7r, nothing happens 
and e/c(7r) = 7T. It is easy to see that in general e^) is a noncrossing 
matching. We refer to the as the Temperley-Lieb operators. (They 
generate an algebra known as the Temperley-Lieb algebra [SI [28], but this 
has no bearing on the present discussion.) 

Theorem 1.8. Define a square matrix H n with rows and columns indexed 
by elements of NC n by 

(H^y = 2n - # {1 < k < 2n : e fe (7r) = tt'} . (5) 

Then /-i n satisfies 

fi n H n = 0. 



Theorems 1.7 and 1.8 seem to be well-known to experts in the field, but 



their precise attribution is unclear to us. As explained in [33] (Sections 



2.2.2 and 3.3.2), Theorem 1.7 follows from the Yang-Baxter equation. A fact 
equivalent to Theorem 1.8 is mentioned without proof in [191 Section 3]. We 
provide a simple proof of this result in Appendix A. 

It is easy to see that the matrix M n = I — ^-H n (where I is the identity 
matrix) is nonnegative and stochastic, i.e., it is a Markov transition matrix, 
associated with yet another Markov chain that has [i n as its stationary vector. 
That is, n n can be computed by solving the vector equation ii n M n = fi n , 
which can be written more explicitly as the linear system 



■in 

2n 

fc=i 

efc(7r') 



X 2 V* (7r e NC n ). (6) 



fc_1 tt'gNC 



The Markov chain (7r m ) m>0 associated with the transition matrix M n has 
the following simple description as a random walk on NC„, known as the 
Temperley-Lieb random walk or Temperley-Lieb stochastic process 

[21J: start with some initial matching tiq] at each step, to obtain 7r m+ i from 
7r m , choose a uniformly random integer k e {1, 2, ... , 2n) (independently of 
all other random choices), and set vr m+1 = efc(7r m ). 

The application of the maps can be represented graphically by associ- 
ating with each a "connection diagram" of the form 
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(a) (b) 

Figure 8: Graphical representation of the application of a sequence of op- 
erators efc, 1 < k < 2n as a "composition of diagrams": (a) the diagrams 
associated with operators e±,ei, e$ are attached to the diagram of the original 
matching; (b) the lines are "pulled" (and any loops are discarded) to arrive 
at the diagram for the transformed matching. 



12 k-1 k k+1 k+2 2n-l 2n 

which will be "composed" with the diagram of the matching tc to which 
is applied by drawing one diagram below the other. (Note that to get the 
correct picture in the case k = 2n one should interpret the diagram as being 
drawn around a cylinder.) By composing a sequence of such diagrams with a 
noncrossing matching diagram one can compute the result of the application 
of the corresponding maps to the matching; see Fig. |8j Using this graphical 
interpretation, it can be seen easily that the stationary distribution n n is 
realized as the distribution of the connectivity pattern of endpoints in an 
infinite composition of e& connection diagrams (where the values of k are 
i.i.d. discrete uniform random variables in {1, ...,2n}), drawn on a semi- 
infinite cylinder. This is illustrated in Fig. |9j 

Yet another interpretation of /-i n is as the ground state eigenvector as- 
sociated with a certain quantum many-body system, the XXZ spin chain. 
The Hamiltonian of this spin chain is an operator acting on the space V = 
(C 2 )® 2 " - , and it has been shown that in the case of "twisted" periodic bound- 
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Figure 9: The connectivity pattern of endpoints in a semi-infinite arrange- 
ment of uniformly random i.i.d. Temperley-Lieb operator diagrams is invari- 
ant under the addition of another operator, hence has \x n as its distribution. 



ary conditions and when a parameter A of the chain, known as the anisotropy 
parameter, is set to the value A = —1/2, the space V will possess a subspace 
U invariant under the action of the Hamiltonian, such that the restriction of 
the Hamiltonian to U coincides (under an appropriate choice of basis) with 
the operator H n defined in ([5]); see [T5J Section 8] and Section 3.2.4]. 

One additional way in which the probability vector n n makes an appear- 
ance is in the study of connectivity patterns associated with a different type 
of loop model known as the fully packed loops (FPLs). An FPL config- 
uration of order n is a subset of the edges of an (n — 1) x (n — 1) square 
lattice [0, n — 1] x [0, n — 1], to which are added 2n of the 4n "boundary" 
edges connecting the square to the rest of the lattice Z, 2 , by starting with 
the edge from (0,0) to (0,-1) and then taking alternating boundary edges 
as one goes around the boundary in a counter-clockwise direction (e.g., one 
would take edges incident to (2, 0), (4, 0), etc.); these extra edges are referred 
to as stubs. The configuration is subject to the condition that any lattice 
vertex in the square is incident to exactly two of the configuration edges; see 



Fig. 10 
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Figure 10: A fully packed loop configuration of order 6 and its associated 
noncrossing matching. 



Let FPL n denote the set of fully packed loop configurations of order n. 
A well-known bijection (22] shows that FPL n is in correspondence with the 
set of alternating sign matrices (ASMs) of order n. These important 
combinatorial objects also have an interpretation as configurations of the 
six-vertex model (a.k.a. square ice) on an n x n lattice with prescribed 
boundary behavior known as the domain wall boundary condition. It 
was conjectured by Mills, Robbins and Rumsey [18] and proved by Zeilberger 
[HT] (see also [3], [TU]) that the number |FPL ra | of ASMs of order n is given 
by the famous sequence of numbers 1, 2, 7, 42, 429, . . ., defined by 

a on T / \ l!4!7!...(3n-2)! 
ASM(n) 



n\(n + 1)! . . . (2n - 1)!' 



As the reader will see below, the function ASM(n) will play a central role 
in our current investigation of connectivity patterns in loop percolation on a 
cylinder. 

Label the 2n stubs around the square [0, n — 1] x [0, n — 1] by the numbers 
1 through 2n, starting with the edge pointing down from (0,0). From the 
definition of fully packed loop configurations we see that any such configu- 
ration induces a connectivity pattern on the stubs, which is a noncrossing 
matching in NC n , in an analogous manner to the way a loop percolation 
configuration on the semi- infinite cylinder does. For 7r e NC n , denote by 
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r\ \ l l n n n ••• r\ 

(a) (b) 
Figure 11: (a) The minimal matching 7r^ in ; (b) the maximal matching 7r™. 



A n (7r) the number of configurations in FPL n whose connectivity pattern is 
equal to n. 

It was first observed by Batchelor, de Gier and Nienhuis [3j that the 
coordinates of the vector related to the enumeration of alternating 

sign matrices. They conjectured the following result, which was later proved 
by Zinn- Justin and Di Francesco [34, 35J. 

Theorem 1.9 (Di Francesco-Zinn- Justin). 1. The numbers /v, (k e NC n ) 
are all fractions of the form a n (7r)/ASM(n) where a n (n) is an integer. 

2. The minimal value of a n (n) is 1 and is attained for n = 7r™ in; the "min- 
imal'^ noncrossing matching consisting of n nested arcs (Fig. 11(a)). 



3. The maximal value of a n (n) is ASM(n — 1) and is attained for = 7r^ a3;; 
the "maximal" noncrossing matching consisting of n nearest-neighbor 



arcs (Fig. 11(b)). 



Shortly after this discovery, Razumov and Stroganov discovered a much 
more precise conjecture (23] about the connection between alternating sign 
matrices and the vector n n ; it turns out that the correct thing to do is to look 
at the ASMs as fully packed loop configurations. Their conjecture, known 
for several years as the Razumov-Stroganov conjecture, was proved in 2010 
by Cantini and Sportiello [5]. 

Theorem 1.10 (Cantini-Sportiello- Razumov-Stroganov theorem). For any 

1 according to a certain partial ordering of noncrossing matchings that is discussed in 
Section [21 
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7r e NC„ we have a n (ir) = A n (n). That is, 

Probabilistically, this means that the vector n n is realized as the distribu- 
tion of the connectivity pattern of a uniformly random FPL configuration of 
order n. Note that this is a finite probability space, whereas the Markov chain 
realization would require a potentially unbounded amount of randomness to 
generate a sample from /j, n . 

Cantini and Sportiello's proof of the conjecture of Razumov and Stroganov 
is highly nontrivial and involves subtle combinatorial and linear algebraic ar- 
guments. Note that while the two (a posteriori equivalent) definitions of [i n 
as the distribution of connectivity patterns in cylindrical loop percolation 
and uniformly random fully packed loop configurations are superficially sim- 
ilar, the two models are quite different. In particular, the cylindrical model 
has an obvious symmetry under rotations of the cylinder by an angle 27r/2n, 
which induces a symmetry under rotation on n n . The fact that the distri- 
bution of the connectivity pattern of random fully packed loops — which are 
defined on a square, not circular, geometry — also has the same symmetry, is 
far from obvious, and its earlier proof by Wieland [29] played an important 
role in Cantini and Sportiello's analysis. 

In this subsection we surveyed several settings in which the probability 
vector ix n appears: as the distribution of the connectivity pattern of cylin- 
drical loop percolation; as the stationary distribution of the Temperley-Lieb 
random walk (or, equivalently, the distribution of the connectivity pattern 
of a semi-infinite arrangement of connection diagrams associated with inde- 
pendent, uniformly random Temperley-Lieb operators); as the ground state 
of the XXZ spin chain; and as the distribution of the connectivity pattern 
of uniformly random fully packed loop configurations of order n. To con- 
clude this discussion, we note that one of the minor results of this paper is 
an additional characterization of the connectivity probabilities (/J. n )weNc n as 



coefficients in a certain natural linear-algebraic expansion; see Theorem 2.12 
in Subsection 12.11 

1.5 Loop percolation on a cylinder: connectivity events 

Having described the origins of the investigations into the random noncross- 
ing matching II* and its distribution fi n , we are ready to discuss the problem 
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of computing explicitly the probabilities of various events. Our main interest 
will be with submatching events of the type {ir < II* n ' ) }, where for tt e NCfc 
and 71 6 NC n (n ^ k) the notation 7r < n means (as defined earlier when n 
is an infinite matching) that tt is a submatching of tt, i.e., that tx\ = 7r . 

|[l,2fc] 

It was observed starting with numerical work of Mitra et al. |19j . Zuber 
[36] and Wilson (unpublished work, cited in |36j) that the probabilities of 
certain events had nice formulas as rational functions in n. For example, one 
has empirically the relations 




3 n 2 + 1 



2 An 2 - 1 ' 

1 97n 6 + 82n 4 - 107n 2 - 792 
8 (An 2 - l) 2 (An 2 -9) 

1 59n 6 + 299n 4 + 866n 2 + 576 
16 (An 2 - l) 2 (4n 2 - 9) 



(n^2), (8) 
(n>2), (9) 



P ( r\ o n (n)N ) 

^123456 ^ * y 



1 214093n 12 -980692n 10 -584436n 8 -1887916n 6 + 13614437i 4 -17432892ri 2 -316353600 . 

= (n > 3 , 

512 (4n 2 -l) 3 (4n 2 -9) 2 (4n 2 -25) 

(10) 

and several other such formulas, which can be discovered by a bit of ex- 
perimentation after programming the linear equations ^ into a computer 
algebra system such as Maple or Mathematical The discovery of additional 
conjectural relations of this type is only difficult insofar as it taxes one's 
patience, programming skill and computational resources, but, as Zuber re- 
marks at the end of [36J, 11 More conjectural expressions have been collected for 
other types of configurations . . . but this seems a gratuitous game in the ab- 
sence of a guiding principle." It therefore appeared sensible to wait for more 
theoretical developments before proceeding with the "gratuitous game." In- 
deed, the first progress (and so far, to our knowledge, the only progress) 
on this front was made by Fonseca and Zinn- Justin [TT] . who, building on 

2 Unfortunately the size of the system is the Catalan number Cat(n) which grows expo- 
nentially with n, making it impractical to compute \i n for values of n much greater than 
n = 9. Zuber }36j managed the computation up to n = 11 by using rotational symmetry 
to reduce the order of the system. 
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a theoretical framework developed earlier by Zinn- Justin and Di Francesco 
[Ml 135] (see also [6]) managed to prove explicit formulas for two classes of 
events. 

Theorem 1.11 (Fonseca-Zinn- Justin [11J). For each k ^ 1, Denote by AC^ 
the "anti- cluster" event that no two of the numbers 1, . . . , k are matched 
under 11* . Denote by B^ 1 the event \7r min < II* [ (the submatching event 



associated with the minimal matching ~n k min ; see Fig. 11). Define a function 



R k (n) 



Then we have 



m 



23-1/ 2 



m 



p(Aci n) ; 



(2j + 1)2)^/2-3 



odd, 



/c even. 



(11) 



P(5 



(n)> 



ASM(A;) i2 fc (fc) 
1 



ASM(ra) 1<ai 



<a 2 <...<a„_ fc 



det 



j + k 
di —j + k 



n~k 



»jj'=l 



In particular, we have P(AC 



2 ; 



5 



2 4n- 



r4j-, which implies equation ^ 



(12) 
(13) 

since 



the event in that identity is the complement of AC 



(n) 



Note that the approach of Fonseca and Zinn- Justin to the relation (12) 



started out, similarly to (13), by expressing the probabilities of these events 



as sums of determinants over a certain family of matrices with binomial 



coefficients. However, in the case of ( 12 ) a family of Pfaffian evaluations due 



to Krattenthaler [15] made it possible to evaluate the sum in closed form. 



As the authors of [TT] point out, both the formulas (12) and (13) have an 
interesting interpretation in terms of the enumeration of certain families of 
totally symmetric self-complementary plane partitions. 



1.6 



Loop percolation on a cylinder: new results and 
the rationality phenomenon 



The identities (|7])-(10) generalize in a straightforward way to a conjecture 
on the form of the dependence on n of probabilities of submatching events. 



21 



Conjecture 1.12 (The rationality phenomenon for submatching events; fi- 
nite n case). IJttq e NC^, the probability of the submatching event Ittq < IT* 1 ^ | 
has the form of a rational function in n, specifically 

p fa < n<->) = - k (" > k ^ ( 14 ) 



ft* i{4n 2_ f)k+ i 



-J 



where Q no (n) is an even polynomial of degree k(k + l) with dyadic rational co- 
efficients. Q no can be computed by polynomial interpolation — see Algorithm B 
in Appendix A. 



In pursuit of an approach that would lead to a proof of Conjecture 1 1 . 12 
one of our goals has been to find a "guiding principle" of the type alluded to 
by Zuber, that would reveal the underlying structure behind the empirical 
phenomena described above. In Section [2] we will prove the following result, 
which can be thought of as one version of such a principle, and is our main 
result. 

Theorem 1.13 (Explicit formula for submatching event probabilities). For 

any tiq e NC&, there exists a multivariate polynomial F no (wi, . . . ,Wk), com- 
putable by an explicit algorithm (Algorithm E in Appendix A), with the fol- 
lowing properties: 

1. F no has integer coefficients. 

2. All the monomials in F 7T0 (wi, . . . , w^) are of the form Oj=i w ^ ^ where 
1 < a\ < . . . < a,}, are integers satisfying ctj < 2j — 1 for 1 ^ j < fcj^ 

3. The probability P ^7r < Tl^^j is given for any n ^ k + 1 by 



p (vro < nf)) = j$^[zi44 ■ ■ ■ z 2 n~ 2 } ( i^Ofc, ...,z k+ 



(15) 



3 Note that the number of different sequences (ai,...,a^) satisfying these conditions 
is Cat(fc), the fcth Catalan number, and indeed the sequence (ai, . . . , au) can be thought 
of as encoding a noncrossing matching in NCfc, a fact that will have a role to play later 
on — see Section [2] 
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where [z™ 1 . . . z™ n ]g(zx, . . . , z n ) denotes the coefficient of the monomial 
z™ 1 . . . z™ n in a polynomial g(zi, . . . , z n ). 



The formula ( 15 ) can be recast in two equivalent forms which some readers 



may find more helpful: first, as a constant term identity 

p fa < n i n) ) = asi^R CT — ( F - ( " 2 ' • • • ' Zk+l) 

x ( z 3 ~ z '^ 1 + z i + Z{Z ^ 



l^i<j^n 1 i-j=l 3 



1 1 7 =1 Z j 



where GT zlr .. )Zn g(zi, . . . , z n ) denotes the constant term of a Laurent poly- 
nomial g(zi, z n ); and second, as a multi-dimensional complex contour 
integral 



P (vro < IlW) = ^g^y f . 4 ( ^ (z 2 , . . . , * fc+1 ) 



X | | ( Z j - Zi)(l + ^ + ZiZj) 2J ._ 1 H 1 — , 

1^2<7=£n 1 1?'=1 / ?'=1 



where the contour for each of the Zj S IS cl circle of arbitrary radius around 0. 

Table [2] lists a few of the simplest submatching events and the multivari- 
ate polynomials associated to them. Note that the first case of the empty 
matching listed in the table (for which the probability of the submatching 
event is 1) corresponds to the known identity 



Zj) 



ASM(n) = \z\z\ . . . z 2 n n - 2 ] Yl & - ^X 1 + z i + z ^ Yl^ 1 + 

(16) 

This identity is already a difficult result. It was proved in |35J by Zinn- Justin 
and Di Francesco, relying conditionally on a conjectural anti-symmetrization 
identity that they discovered which itself was proved shortly afterwards by 
Zeilberger [32]. Their proof also relies on a highly nontrivial Pfaffian evalu- 
ation due to Andrews [1] (see also [21 US])- The following conjecture can be 



thought of as a natural generalization or "deformation" of (16) 
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k 


7T 


W\...W2k 


n 
u 


empty 


1 
1 


1 




1 


1 2 




2 




1 




12 3 4 




2 








12 3 4 





7T 



^ ^ ^ 

1 2 3 4 5 6 



1 2 3 4 5 6 



1 2 3 4 5 6 



1 2 3 4 5 6 




UJ i...W2fc 



W 2 - WxW2 



W 1 



WiW 2 



Table 2: The polynomials F no associated with some submatching events, 
after factoring out the product wi . . . w 2 k which always divides F no . (The 
first case of the "empty" matching refers to the trivial matching of order 
k = 0, in which case the submatching event has probability 1.) 



Conjecture 1.14 (The algebraic rationality phenomenon). Let k ^ 1, and 

let 1 ^ ai < . . . < a*; be integers satisfying aj < 2j — 1 for all j . There exists 
a rational function of the form R(n) = P(n)/ nj=i(4^ 2 — j 2 ) k+l ^ where 
P{n) is an even polynomial of degree at most k{k + l)/2 with dyadic rational 
coefficients, such that for all n ^ k + 1 we have that 

/ k n \ 

■ ■ ■ 4 n - 2 ] n ^ n & - + ^ + n (! + ^ 

= ASM(n)i2(n). (17) 
It is unclear whether the specific assumption about the sequence ai, . . . , 
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that enters the form of the monomial Ylj=i z ^+\' 3 ls necessary to imply the 



rationality of R(n) in (17). The main assumption, which may be sufficient 



to imply the result, is that we are looking at a Taylor coefficient that is "a 



fixed distance away" from the coefficient in (16) 



An immediate consequence of Theorem 1.13 is that Conjecture 1.14 



es- 



sentially implies Conjecture |1. 12 although a small gap remains regarding the 
validity of (14) in the case n = k. 



Theorem 1.15. Conjecture \l . 14\ implies a weaker version of Conjecture 1.12 
in which ( 14 ) is only claimed to hold for n 5= k + 1 . 



From the above discussion we see that, while we have been unable to get a 
complete understanding of the probabilities of submatching events, we have 
reduced the problem to the essentially algebraic question of understanding 
the form of the dependence of n of the polynomial coefficients appearing on 



the left-hand side of (17). Moreover, the algorithm for finding the polynomial 



F 1TQ associated with a submatching event, which will be explained in Section [2j 
stems from a fairly detailed theoretical understanding of the model, and 
therefore already helps eliminate much of the mystery surrounding identities 
such as 



It should be noted as well that the relations (16) and (17) belong to a 



large family of identities known as constant term identities. The study 
of such identities became popular following the discovery by Dyson of the 
identity 



CT 2 




{a\ + . . . + a r 
ai! . . . a n \ 



(ai,.. .,a n > 0), 



(18) 

which became known as the Dyson conjecture [7]. (Dyson's conjecture was 
proved by Gunson [T4] and Wilson [30] , and a particularly simple proof was 
later found by Good [12].) Research on such identities has been an active 
area that involves a mixture of techniques from algebraic combinatorics and 
the theory of special functions. In particular, Sills and Zeilberger [25] proved 
a deformation of Dyson's identity in which a coefficient near the constant 
term is shown to equal the multinomial coefficient on the right-hand side of 
( 18 ) times a rational function in the exponents a\, . . . , a n ; this is quite similar 



in spirit to the claim of Conjecture |1.14 



To conclude this section, we show how the relation (|8j) can be derived in a 
simple manner from the results of Fonseca and Zinn- Justin, and in addition 
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derive another identity concerning a finite connectivity event which is not a 
sub matching event. 

Theorem 1.16. The identity <^ holds for n ^ 2^ and we have the addi- 
tional identity 

15 (n 2 -4)(9n 4 + 38n 2 ~63) 
" 16 ' (4n 2 - l) 2 (4n 2 - 9) 1 " J 1 J 

fwe use i/je notation for submatching events for convenience, but note that 

in) n (n> n (n) 

this refers to the event that II* matches the two pairs 1 <-^-> 2 and 4 <-^> 5 / 



Proof. By (12), we know that 

FrAr (n),_33 (n 2 -l)(n 2 ~4)(n 2 -9) 

F(A ° 4 j " ¥ ' (4n 2 -l) 2 (4n 2 -9) • (20) 

On the other hand, by the inclusion-exclusion principle, we have that 

p(Acl n) ) = i - p e yy) ~ p ( n * n) e O ) ~ p ( n * rt) e 




(using ([7]) and the rotation-invariance of U^). Substituting the probability 
from Q and solving for P (lli n) e f \ \ gives (pBI). For the second 
identity (19), perform a similar inclusion-exclusion computation for the event 
AC5 , whose probability is given according to ( Jl2| by 

PrAr (^ = H (n 2 -l)(n 2 -4)(n 2 -9) 
1 5 j 16 (4n 2 -l) 2 (4n 2 -9) ' 

The details of the computation are easy and left to the reader. □ 



4 This proves part of Conjecture 9 in Zuber's paper 
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1.7 Consequences for loop percolation on a half-plane 

Let us return to the original setting of loop percolation on a half-plane dis- 
cussed earlier. The following result allows us to deduce exact results on 
probabilities of local connectivity events in the half-plane from correspond- 
ing results in the cylindrical model. 

Lemma 1.17. Let A a Z x Z be a finite set, and let E n denote the finite 
connectivity events 



e.- n {;^4 

Co Me A V J 



(j,k)eA 

associated with A (which are defined for large enough n). Then we have 
F(E n ) - F(E) asn^cc, where E = C\ Utk)eA {j " k }- 

Proof. Couple the cylindrical and half-planar loop percolation by using the 
same random bits to select the edge configuration incident to the vertices in 
the region V(LP n ) (defined in (|3|). With this coupling, it is easy to see that 
the symmetric difference E n AE of the events E n and E is contained in the 
event B n that one of the half-planar loop percolation paths starting at the 
point (j, — j) for some j belonging to one of the pairs (j, k) e A reaches the 
complement Z 2 \V(LP n ). Since the loop percolation paths are almost surely 
finite, we have \im n ^ X) F(B n ) = P (p|™ =1 -B n ) = 0, and therefore we get that 
\F(E n ) - F(E)\ «S F(E n AE) ^ F(B n ) -> as n -> oo. □ 

As immediate corollaries from the lemma we get the following facts: first, 



Conjecture 1.12 (or the weaker version of it mentioned in Theorem 1.15) 



implies Conjecture 1.3 Second, Theorem 1.4 follows as a limiting case of 



([7]) and (|8p. Third, we have the following explicit formulas for events in the 
half-plane model corresponding to (12) and (19). 



Theorem 1.18. For k ^ 1, let ACk denote the anti-cluster event that no 
two of the numbers 1, . . . , k are matched under II*. We have the formulas 

F(AC k ' 



2lfc/2jW2+ljASM(A;)i^ fc (A;) , 
where Rk(-) is defined in (11) (see Table^, and 



135 
1024' 
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k 


i 


2 


3 


4 


5 


6 


7 


p(n*eAcy 


5 


1 


33 


11 


2431 


85 


126293 


8 


4 


512 


1024 


221 


220 


2 35 



Table 3: Probabilities of the anti-cluster event ACu for k 



The theory of wheel polynomials and the 
qKZ equation 



Our goal in this section is to prove Theorem 1.13 The proof will be based 
on an extension of an algebraic theory that was developed in a recent series 
of papers [TTI I3"3"j [Ml ES], where it is shown that a tool from the statisti- 
cal physics literature known as the quantum Knizhnik-Zamolodchikov 
equation (or qKZ equation) can be applied to the study of the connectivity 
pattern of cylindrical loop percolation. 

Mathematically, the qKZ equation as applied to the present setting re- 
duces to the analysis of a vector space of multivariate polynomials satisfying 
a condition known as the wheel condition. We call such polynomials wheel 
polynomials. In Subsections 2.1 - |2.3 we will survey the known results re- 
garding the theory of wheel polynomials that are needed for our purposes. 



In Subsection 2.4 we present a new result (Theorem 2.24) on expansions of 
families of wheel polynomials associated with submatching events, that is 
another main result of the paper. In Subsection 2.5 we show how to derive 
Theorem 11.131 from Theorem 12.241 



2.1 Background 

Let NC n denote as before the set of noncrossing matchings of 1, . . . ,2n. It 
is well-known that elements of NC n are in canonical bijection with the set 
of Dyck paths of length 2n, and with the set of Young diagrams contained 
in the staircase shape (n — 1, n — 2, . . . , 1). These bijections will play an 



important role, and are illustrated in Fig. 12 ; see [TH Section 2] for a detailed 
explanation. 

We adopt the following notation related to these bijections. If 7r e NC n , 
let Tik be 1 if n(k) > k ( ll k is matched to the right") or —1 if ir(k) < k ( ll k 
is matched to the left"). The vector (iri, . . . , 7r 2n ) is the encoding of 7r as the 
sequence of steps in the associated Dyck path. Let 7r + denote the sequence 
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123456789 10 123456789 10 

Figure 12: The noncrossing matching from Fig. E[a) represented as a Dyck 
path and as a Young diagram. 



(oi, . . . , a n ) of positions where 7Tfc = 1 (the positions of increase of the Dyck 
path). Encoding ir in this way maps NC n bijectively onto the set 

NC+ = {a = (ai, . . . ,a n ) : 1 < a x < . . . < a n < 2n - 1 

and dj < 2 j - 1 for all j}. (21) 

The Young diagram associated with a noncrossing matching 7r e NC n is 
denoted X n . For n, o e NC n , denote <r /" 7r if A CT is obtained from by the 
addition of a single box, and it < o if A,r is contained in A CT . (This partial 
order on NC n is precisely the order with respect to which the matchings 7r^ in 



and vr^ ax defined in Subsection 1.4 are minimal and maximal, respectively.) 



Denote n /*,• a if X a is obtained from A T by adding a box in a position that, 



in the coordinate system of Fig. 12, lies vertically above the positions j and 
j + 1 on the horizontal axis. 

Definition 2.1 (Wheel polynomials). Let q = e 2m ^. A wheel polyno- 
mial of order n is a polynomial p(z) = p(z\, . . . , Z2 n ) having the following 
properties: 

1. p is a homogeneous polynomial of total degree n(n — 1); 

2. p satisfies the wheel condition 



p( Zl ,...,z 2n )\ =0 (Ui<j<K2n). (22) 

Denote by W n [z] the vector space of wheel polynomials of order n. 
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It is easy to see that the condition (22) depends only on the cyclical order- 
ing of the variables Z\, . . . , Z2 n (in other words, the space lf„[z] is invariant 
under the rotation action p(zi, . . . , Z2 n ) >-*■ p(z 2 , . . . , Zi ni 21)); hence the name 
"wheel polynomials." 

The algebraic properties of wheel polynomials are quite elegant and of 
direct relevance to our study of connectivity patterns in loop percolation. 
Below we survey some of the known theory. 

Example 2.2. Let p(z) = ]li^<K« ~ ^ z i) Yl n +i^i<j^2n (<l z i ~ 
It is easy to check that p is a wheel polynomial. Indeed, it is homogeneous of 
the correct degree, and if 1 < % < j < k < 2n and we make the substitution 
Zk = q 2 zj = q 4 Zi, then, if j < n there will be a zero factor in the first product 
Q 1<i<j<n (qzi — q~ l Zj) 1 and similarly if j > n then there will be a zero factor 
in the second product. In either case the wheel condition (22) is satisfied. 

Two types of pointwise evaluations of a wheel polynomial will be of partic- 
ular interest: at the point z = (1, 1, . . . , 1) and at a point z = (g _7ri , . . . , q~^ 2n ) 
associated with the steps of a Dyck path, where 7r e NC n . As a shorthand, 
we denote p(l) = p(l, . . . , 1) and p(ir) = p(q~ ni , ■ ■ ■ , q^ 2n ) and refer to these 
values as the 1-evaluation and 7r-evaluation of p, respectively. 

Theorem 2.3. A polynomial p e VT^fz] is determined uniquely by its values 
p(tt) as 7r ranges over the noncrossing matchings in NC n . Equivalently, the 
Ti -evaluation linear junctionals (ev 7r ) 7r6 Nc n defined by ev^p) = p(n) span the 
dual vector space (W n [z])*. 

Proof. See [TUl Appendix C]. □ 

In particular, it follows that dimM/^fz] < |NC n | = Cat(n). We shall soon 
see that this is an equality, that is, the evaluation functionals (ev 7r ) 7reNCn 
are in fact a basis of (W n [z])*. The proof of this fact involves a remarkable 
explicit construction of a family of wheel polynomials that will turn out to 
be a basis of Wnfz] dual to (ev 7r ) 7re Nc„, and will play a central role in our 
analysis. 

For a multivariate polynomial p e C[z\, . . . , z m ] and 1 < j < m, define 
the divided difference operator dj : C[zi, . . . , z m ] — > C[zi, . . . , z m ] by 



(d j p)(z 1 ,...,z m ) 



p(Zi, . . . , Zj + i, Zj, . . . , Z2n) — P[Zi, . . . ,Zj, Zj + i, . . . , Z 2n ) 



Zj+i Zj 

The next two lemmas are easy to check; see [531 Sec. 4.2.1] for a related 
discussion. 
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Lemma 2.4. If p e W n [z] and 1 < j < 2n — 1, £/ien (gzj — g 1 Zj- + x)5_ ? -p 
a/so m lf n [z]. 

Lemma 2.5. If for some 1 < j < 2n — 1, it e NC n satisfies j * — » j + 1, 
£/jen t/ie set of preimages ej 1 (7r) o/7T under ej (t/ie Temperley-Lieb operator 
defined in (jil) ) consists of: 



1. a noncrossing matching n e NC n satisfying it / 'j it, if such a matching 
exists (i.e., if the associated Young diagram is still contained in the 
staircase shape (n — 1, . . . , 1)); and 

2. a set of noncrossing matchings a e NC„ satisfying a < it. 

The last two lemmas make it possible to construct a family of wheel poly- 
nomials indexed by noncrossing matchings in NC n (or, equivalently, by Young 
diagrams contained in the staircase shape (n— 1, . . . , 1)) recursively. We start 
with an explicit polynomial known to be in H^z] — a scalar multiple of the 
polynomial from Example |2.2 above — which we associate with the minimal 
matching 7r™ in (which corresponds to the empty Young diagram). We then 
define for each noncrossing matching 7r a polynomial obtained from the poly- 
nomials of matchings preceding it in the order < using linear combinations 
and the operation of Lemma 2A The precise definition is as follows. 

Definition 2.6 (qKZ basis). The qKZ polynomials are a family of poly- 
nomials (^TOvreNCn defined using the following recursion on Young diagrams: 

%njz) = (-3)-© Yl {qzi-q- 1 ^) f] (g* - q-'zj) , (23) 
* w (z) = (qzj - q~ 1 z j+1 ) dj^ a - Yi ^ ifa/jTT. (24) 

Since for a given noncrossing matching n, the choices of a and j such that 
a / 'j 7T are in general not unique, it is not clear that the above definition 
makes sense. The fact that it does is a nontrivial statement, which we include 
as part of the next result. 

Theorem 2.7. The qKZ polynomials satisfy: 

1. ^ n is well-defined, i.e., the result of the recursive computation does not 
depend on the order in which boxes are added to the Young diagram. 
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2. tt ff e W n [z]. 



3. Let p(ix) denote the rotation operator acting on tt, defined by (p(jr))(k) = 
ix{k + 1) for 1 < k < 2n and (p(ir))(2n) = vr(l). Then we have 

^ P ( 7T )(zi,...,z 2n ) = ^ 7V (z 2 ,...,z 2n ,z 1 ). (25) 
Proof. See [331 Section 4.2]. □ 



As explained in [331 Sections 4.1-4.2], equations (23)-(25) together are 
equivalent to a different system of equations which forms (a special case of) 
the qKZ equation. A different way of solving the same system, which we will 
not use here, is described in [9]. 

If 7i e NC n and 1 ^ j ^ 2n — 1 is a number such that j <-^-> j + 1, 
denote by ftj the noncrossing matching in NC n _i obtained by deleting the 
arc connecting j and j + 1 from the diagram of the matching and relabelling 
the remaining elements. This operation makes it possible to perform many 
computations recursively. The next result provides an important example. 

Theorem 2.8. If j j + le NC n then for any a e NC n we have 

, , f3 n - 1 ^.(0- 7 -) if 7 "7 + 1, , , 

V v {a) = { A 3> J J (26) 

1 otherwise. 

Proof. See [331 Section 4.2]. □ 
Corollary 2.9. i. For a// tt, a e NC n , we /iawe 

*» = <W = {* ^ = CT ' (27) 
1 otherwise. 

2. The qKZ polynomials {^t^-kenc^, are a basis for W n [z] . 

3. The n-evaluation linear functionals (ev 7r ) 7re Nc„ are the basis o/(W n [z])* 

^. dimW n [z] = Cat(n). 



Proof. Part 1 follows by induction from (26), and parts 2-4 are immediate 



from part 1. □ 
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Example 2.10 ( [331 Section 4.3.1]). Let p(z) = s\(zi, ■ ■ ■ , z 2n ), the Schur 
polynomial associated with the Young diagram A = (n— 1, n— 1, . . . , 2, 2, 1, 1). 
That is, we have explicitly 

P(z)=( EI det^)^ (28) 

\l!g2<j!£2n / (Kj^3n-2, j=0,l (mod 3) 

Then p is homogeneous of degree n(n — 1), and the matrix whose determinant 



appears in (28) has the property that if we make the substitution 



q Zj = q Zi for some l^i<j<k^ 2n, the columns with index i,j, k 
of the matrix become linearly dependent (take the linear combination with 
coefficients l,q 2 , q A ), causing p to be 0. Thus, p is a wheel polynomial. 

Since p is also a symmetric polynomial, it satisfies p(vr) = p(7T ; ) for any 



7r, 7r' g NC n . By Theorem 2.3 it follows that p is the unique symmetric wheel 
polynomial up to scalar multiplication. Note that this Schur function is also 
known to be the partition function of the six-vertex model with domain wall 
boundary condition at the "combinatorial point," and played an important 
role in the study of the enumeration of alternating sign matrices; see [201 [27] 
and [331 Section 2.5.6]. 

Denote ip n = ^(1) and xj) n = (t/v^sncv The vector tp n is related to the 
connectivity pattern in cylindrical loop percolation, as explained in the 
following result. 

Theorem 2.11. The vector tj) n is a solution of the linear system ([61. Conse- 
quently it is a scalar multiple of the probability distribution fi n q/TI* . That 
is, for any ir e NC n we have ipn = a n^n, where a n = Xi7reNC„ 

Proof. See [331 Sec. 4.3]. □ 



We will see later (see Theorem 2.22 below) that the normalization con 



stant a n is equal to ASM(n). Assuming this fact temporarily, Theorem 2.11 
implies the following interesting and purely algebraic characterization of the 
numbers = ¥(11^ = n), which does not seem to have been noted before. 

Theorem 2.12. For any wheel polynomial p e W n [z], we have that p{l) = 
ASM(n) Xi7reNC„ A i i-P( 7r )- In other words, the 1- evaluation functional ev x has 
the expansion 

evx = ASM(n) ev vr 

7TGNC n 
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as a linear combination of the ir-evaluation Junctionals. 

Proof. If p = XivreNCn ^t^t gives the expansion of p in terms of the qKZ basis 
polynomials, then 

p(l)= J] 6A(1)= 2 Mv = ASM(n) J] 

TTSNCn TTSNCn 7rSNC„ 

On the other hand, taking the ^-evaluation of the expansion p = 2o-eNC„ ^>^a 
easily implies using (27) that b n = p(n). □ 

Next, we construct a second important family of wheel polynomials de- 
fined using contour integration. 

Definition 2.13 (Integral wheel polynomials). Denote 

A n = {a = (ai, . . . , a n ) : 1 < ai < . . . < a n < 2n — 1 

and a,j < 2j — 1 for all j } . 

(These sequences generalize the steps of increase of Dyck paths — compare 
with (21 ).) For a sequence a = (a 1; . . . , a n ) e A n , denote 

$ a (zi, . . . ,Z 2 n) = Y\ ^ ~ ^Zfc) 

Y] (wk - Wj)(qvjj - q~ lr w k ) 




where the contour of integration for each variable Wj surrounds in a coun- 
terclockwise direction the singularities at z fc (1 ^ k < aj), but does not 
surround the singularities at q~ 2 Zk (% + 1 < k < 2n). 

We will be mostly interested in <3> a for a e NC^ — the set of such polyno- 
mials will turn out to be a basis for W n [z] — but the case when a is in the 
larger set A n will play a useful intermediate role. 

The above definition of $ a is somewhat difficult to work with. In par- 
ticular, it is not even obvious that the $ a 's are polynomials. However, the 
contour integrals can be evaluated using the residue formula, which gives a 
more concrete representation of the $ a 's that makes it possible to show that 
they are in fact wheel polynomials and prove additional properties. 
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Lemma 2.14. 1. The functions $ a can be written explicitly as follows. 
$ a (z) = (-l)© Yl (q Zj -q-'z k ) 

l^j<k^2n 



2 



(m i , . . . , m Tl ) 

1 )Tlj Oj , 771 j ^ m/. 



n 

3=1 



2?i 



]~[ (2 m , - 2*) 1 [ (qz mj - q l z k ) 



\ k ^ a,j 



Zk) 



2 

Af — (mi , . . . , mn ) 



j\inv(M) 



1 ^ j < fe ^ 2n 
j £ M or j — rrip, k sr ; a p 



j — l 1 ^ fe ^ a_j 

fe $ iW ox fe > m-i 



where for a sequence M = (mi, . . . ,m n ), inv(M) denotes the number 
of inversions of M , given by inv(M) = < z < j < n : > Jn^}. 



$ a a u>/jee/ polynomial of order n. 
Proof. See [53 Section 3]. 



□ 



For each sequence a = (ai,. . . , a n ) e „4 n , let (C aiCr ) - e Nc n be coefficients 
such that 

o-eNC„ 

For 7r e NC n denote C %>(r = C a)Cr with a = n + . The coefficients (C 7rj(T ) - e Nc„ 
have the property that 

o-eNC„ 

Lemma 2.15. We have the relations 



®*(o) = C gL>a (as A n ), 
$M=C V><T (tteNC„). 
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Proof. This follows immediately from (27) by considering the 7r-evaluation 
of both sides of □ 

The next result is a recurrence relation that makes it possible to recur- 



sively compute the quantities $ a (o") = C a ,cr, analogously to Theorem 2.8 



Theorem 2.16. If a s NC ra has a "little arc" j <-^-> j + 1 and a e A n , let p 

be the number of times j appears in a ; so that we can write 

p times 

a = (ai, . . . , dk,j, ■ ■ ■ , j, Ofe+p+i, • • • , a n ) 
where < j < ctfc+p+i- In the case when p > 0, denote 

p— 1 times 



(Oi, 



• a>k,J, 



,j,Cbk+p+l 2 



, . . . , u< n 



TTien we /tave 
where 



3 n 1 x(p)®z j {v j ) ifp>0 (i.e.,ifje{a 1 ,...,a n }), 
ifp = 0, 



(32) 



(33) 



x(p) 



q p - q- p 
q - q- 1 



ifp = 0(mod3), 

1 if p = 1 (mod 3), 
— 1 if p = 2 (mod 3). 



(34) 



Proof. See [351 Appendix A]. 



□ 



Note that if a = vr + e NC+ then a,- is in A n but not necessarily in NC^, 
which is the reason why it was necessary to consider the family of polynomials 
indexed by the larger set A n even though we are mainly interested in the 
polynomials $„■. 

Like (26), the recurrence relation (33) can be solved to give an explicit 
formula for C aa . 

Theorem 2.17. For a = (ai, . . . ,a n ) e NC^ and o e NC n; the coefficient 
Ca,o- is given by: 

C a>(T = n x(Pi.fc(a)), (35) 

jr'<fc, j*^k 

where we denote Pj,k(&) = #{ m : j ^ a m < A;} — — j — 1). 
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Proof. See (33 Appendix A]. □ 

Theorem 2.18. The matrix (C 7r . -) 7riCrG Nc„ is triangular with respect to the 
partial order < on link patterns, and has Is on the diagonal (i.e., = 1 
for all 7r e NC ny ), hence it is invertible. 

Proof. See HH Sec. 3.3.2]. □ 

Corollary 2.19. ( < & 7r ) 7rG NC, l form a basis for W n [z]. 

Denote the coefficients of the inverse matrix of (C 7I - iCr ) 7ri(7e Nc„ by Cy.o-j or 
C^a if a = cr + . In other words, the coefficients C n ^ a are defined by the 
relation 

Y[r = V C $ 

creNCn 

We have defined two square matrices C n = (C 7ri0 -) vr)Cre Nc, l an d C n = C^ 1 = 
(Cn,a)-K,ae'Nc n , which are the change of basis matrices transitioning between 
the two wheel polynomial bases we defined. Note that C n%(T , like C n>a , are 
integers, but there does not seem to be a simple formula for them. The 
matrices C n and C n for n = 2,3,4 are shown in Table [4} As we shall see, 
these matrices will enable us to translate certain sums of the qKZ polynomials 
(which are associated with certain connectivity events in loop percolation) to 
linear combinations of the ^'s. To get back to numerical results about the 
numbers ?/v = ^(1) (which are related to the probabilities fi n according to 
Theorem 2.11), we will want to consider also the 1-evaluations of the ^v's. 
Denote a = $ a (l), 0vr = ^(l)- The coefficients 4> w are related by 

creNC n 

^ = Zj ^afa. 
creNC„ 

Lemma 2.20 ([35J, Section 3.4). a has the following contour integral for- 
mula: 



1 du 

Y\ (tt fc - Uj)(l + u k + UjU k ) Y[-aj- ■ Oo - ^"' ( 36 ) 

l=£j<fc=£n 3=1 3 3=1 



where the contours are circles of arbitrary radius around 0. Equivalently, 
can be written as a constant term 



\l^j<k^n 3=1 U j J 
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Table 4: The matrices C n , C n for n = 2,3,4 (the ordering of the rows and 
columns corresponds to decreasing lexicographic order of the Young diagrams 
associated with noncrossing matchings). Note that C n can have entries larger 
in magnitude than 1; the smallest value of n for which this happens is n = 7. 



or as a Taylor coefficient 



[ul 1 1 ... <" X ] I Y\ ( u k - + u k + UjU h ) J . 



Proof. In the integral formula ( 29 ) set Z\ 



■ ■ ■ = = 1 and then make the 
substitution Wj = 1 ^_ 9 gM " J ■ This gives (36) after a short computation. □ 



2.2 Product-type expansions and the sum rule 

We are ready to start applying the theory developed above to the prob- 
lem of computing probabilities of connectivity events for the random non- 
crossing matching II* n \ Our methods generalize those of Zinn- Justin and 
Di- Francesco [35], so let us recall briefly one of their beautiful discoveries. 
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Define a set B n <= NC+ consisting of 2 n 1 sequences, given by 

B n = {a = (oi, . . . , O e An : a% = 1, a 3 - e {2j - 2, 2j - 1}, (2 < j < n)}. 
For each a = (ai, . . . , a n ) e B n define 

£(a) = {tt e NC n : for all 1 < j < n, 7r 2i _i = 1 iff Oj = 2j - 1}. 



Note that NC„ = |J £( 



a . 



aeB n 



Theorem 2.21. For any a e B n , we have the expansion 

7re£(a) 



(38) 



Theorem 2.21 and its elegant proof served as the inspiration for one of our 



main results (Theorem 2.24 in Subsection 2.4). As an aid in motivating and 
understanding our use of Zinn- Justin and Di Francesco's proof technique in 
a more complicated setting, we include their original proof, taken from [35| 
Section 3], below. 



Proof of Theorem 2.21, By Theorem 2.3, it is enough to prove that for any 



a e NC„ we have 



s 

7re£(a) 



(39) 



The proof is by induction on n. Denote the right-hand side of (39) by 9 a (a 



Let k <^~* k + 1 be a little arc of a. Let j be such that k e {2j — 2, 2j — 1}. 
By the definition of B n , either k = aj is the only occurrence of k in a, or k 
does not occur in a (in which case aj is the number in {2j — 2, 2j — 1} of the 
opposite parity from k). For the left-hand side of (39) we have by (33) that 



jn— 1 



$a fc (o-fc) ifa,- 



k. 



if aj k. 



From the definition of B n it is easy to see that when k = aj, slj- e B n -\. So, 
to complete the inductive proof it is enough to show that 6 a (a) satisfies the 
same recurrence. First, assume that aj k. If k = 2j — 1 (so aj = 2j — 2) 
then any n e £(a) will have the property that 71^ = —1, that is, n matches k 
with a number to the left of k. In particular 7r does not contain the little arc 
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(k, k + 1), and therefore ^(cr) = 0. This shows that # a (cx) = 0. Alternatively, 
if k = 2j — 2,dj = 2j — 1, then any 7r e /3(a) will have the property that ir 
matches aj = k + 1 with a number to the right of k + 1, and again the little 
arc (fc, fc + 1) cannot be in tc and ^(cr) = 0, so the relation 9 a (a) = also 
holds. 

Finally, in the case aj = k, the noncrossing matchings n e /2(a) can 
be divided into those for which (k, k + 1) is not a little arc of it, which 
satisfy ^ n (cr) = 0, and those for which (k,k + 1) is a little arc of it. The 
ones in the latter class satisfy ^(cr) = S™ -1 ^^^), and furthermore the 
correspondence n i-> n k maps them bijectively onto £(a fc ). Thus we get the 
relation 

e a (a) = 3 n - 1 e^ k (a k ), 

which is the same recurrence as the one for $ a ( cr ) an d therefore completes 
the proof. □ 



Theorem |2.21 has several immediate corollaries that are directly relevant 



to the problem of deriving formulas for probabilities of connectivity events. 



First, summing (38) over the sequences aeB n gives that 



2 ».-2*- 

7TGNCn aeBn 

Taking the 1-evaluation of both sides then gives the relation 

7reNC„ aeB„ 



and the right-hand side can then be expanded as a constant term using (37), 
which gives that 

J] = CT Uli ... )Un j Yl ( u k ~ + u k + UjU k ) Y[ 1 tj-2 ) • ( 40 ) 

Note that the key point that makes such an expansion possible is the fact 
that B n decomposes as a Cartesian product, which manifests itself as the 
product 0^=2(1 + u j)/ u2 /^ 2 i n the Laurent polynomial. It is precisely this 
product structure that our more general result seeks to emulate. Finally, 



applying the identity (16) gives the following result. 



Theorem 2.22. EtteNC Vv = ASM(n). 



40 



With this result, we are finally able to relate results about the vector t/> n 
to statements about its normalized version, the probability vector n n . As 
an example (also taken from [35]), letting a = (1, 3, 5, . . . , In — 1), note that 
£(a) is the singleton set {vr" ax }, so, applying (38) in this case, then taking 



the 1-evaluation and using (|37|) as before, gives the identity 

Ip-jj-n = CT,, 
' max 

This constant term can be evaluated iteratively as 



- Ul ,...,un I Yl ( Uk ~ + Uk + u i Uk ^ fl ) 



CT ttl) ... )tln _ 1 CT UB ( Yl K - + u k + UjU k ) Y\ ^2 ) 

( n 1 + U 

= CT Ul) ... iU „_ 1 j Yl (U k - Uj)(l + U k + UjU k )Yl 2 j-2 

3=1 U J 



K l^j<k^n-1 



ASM(n - 1 : 



(use (16), noting that on the right-hand side of that identity, starting the 
product Yljfi + z j) at j = 1 instead of j = 2 does not affect the constant 



term), confirming part 3 of Theorem 1.9 



2.3 ^-nested matchings 

If 7r e NC n , denote by JT the noncrossing matching of order n + 1 obtained 
from 7r by relabeling the elements of [1, 2n] as 2, . . . , 2n + 1 and adding a 
"large arc" connecting 1 and In + 2. We call this the nesting of n or u tt- 
nested." For an integer p ^ let JT P denote the p-nesting of tt, defined as 
the result of p successive nesting operations applied to tt. Similarly, for the 
associated sequence a = n + e NC^, denote IT = (^T) + and If p = ("7T P ) + . 
If a = (ai, . . . , a n ), it is easy to see that this can be written more explicitly 
as 

•£T P = (1, . . . ,p,p + ai, . . . ,p + a n ). 
The following result is due to Fonseca and Zinn- Justin [TT] . 
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Theorem 2.23. Let p ^ 0. Under the p-nesting operation, the expansions 

7reNC n 

2 (aENC n ), 



\1> 

* 7T 



CTSNCn 



tnrn into 



TTSNCn 



o-eNC 



Equivalently, we have the recursive relations 



otherwise, 



otherwise. 



(41) 
(42) 



Proof. The relation (42) is an immediate consequence of (41). To prove 



(41), note that the fact that C^p^ = if /i is not of the form a p for 



some a e NC n follows from the triangularity of the basis change matrix with 
respect to the order <. The relation C-^ P ^v = C CTi7r is easy to prove by 
induction from (33). □ 



2.4 General submatching events 

Fix k ^ and tt e NC fc . For n ^ k + 1 define 

£ n (ir ) = {vr g NC n : tt(j) = 1 + tt (j - 1) for 2 < j < A; + 1} . 



In the terminology of Subsections 1 1 .3| and 1 1 . 5[ this set of matchings is related 
to the submatching event associated with the submatching 7r , except that we 
require the submatching to occur on the sites in the interval [2, k + 1] instead 
of [l,k]. Surprisingly, this turns out to be the correct thing to do — see below. 
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Next, for any n ^ k+1 and b = (bk+2, ■ ■ ■ , b n ) satisfying bj e {2j— 2, 2 j — 1 } 
for all j, define 

£ n (ir , b) = iir 6 £ n (7r ) : 7r 2i _i = 1 iff = 2j - 1 (k + 2 < j < n)|. 



Note that 



^n(TTo) 



|J £:„(7r ,b). 



b=(6&+2,— ,&n) 

Vj 6 J -e{2i-2,2j-l} 



The following result will easily imply Theorem |1.13[ and is one of the 
main results of this paper. 

Theorem 2.24. For any tc g NC&, n ^ k + 1 and vector b = (&fc+2, • • • , & n ) 
satisfying bj e {2j — 2, 2j — 1} /or a// j, we /jcrae 



/j ^ ^7ro,a^(l,l+a 1 ,...,l+a fc ,& fe+ 2,. ..,&«)■ 

7r££ n (7r ,b) a=(ai,...,a fc )sNC+ 

It follows that for any ttq e NQ. and ra ^ + 1, 



(43) 



7re£„(7ro) 



l,l+oi,...,l+Ofc,6fc + 2,...,6. 



>„)■ (44) 



b k+2 ,...,b n a=(ai,...,a fc )eNC+ 
Vj b j6 {2i-2,2i-l} 



Proof. By (42), the right-hand side of (43) is equal to 



C TO ,a < ^ > (a 1 ,...,a fe+1 ,6 fc+2 ,...,fe„) • 



a=(ai,...,a fe+1 )eNC+ +1 



So, instead of proving (43) we will show that 



yi ^ C vro ,a C &(ai,...,a fc+1 ,6 fe+ 2,...,fen)- (45) 

7T6£ n (7r ,b) a=(ai a i+ i)eNC+ +1 

We will prove by induction on n the following slightly stronger claim: if we 
have an expansion of the form 



2 C a$a 

ae.4 fc+ i 
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(46) 



for some constants c a (a e Ak+i), then for any vector b = (b k+2 , ■ ■ ■ ,b n ) 
satisfying bj e {2j — 2,2j — 1} as in the theorem, the equality 

Yj **= S C a^(a,b) (47) 
7rGf„(7ro,b) ae_4 fc+1 



holds. (Note that this claim implies (45) by taking c a = C-^a for a e NC+ 



and c a = otherwise, but allowing the sum on the right-hand side of (46) 



to range over the bigger set A k +i helps make the induction work.) Denote 



the left-hand side of (47) by X no>h and the right-hand side by Y ch . We now 
prove that = Y c ^. 

The induction base: here k = 0, n = 1, 7Tq = (the empty matching), 
Ak+i = A\ = {1} and the expansion necessarily takes the form 



Also £ n (7r , b) = £((1, b 2 , ■ ■ ■ , 6 n )), so the claim reduces to the relation (38) 
from Theorem 12.211 



The inductive step: First, if n = k + 1 then the claim (47) is the same 



as the assumption (46), since b is a trivial vector of length and £n(ffo, b) is 
the singleton set {'vf^}, so there is nothing to prove. From now on assume 
n > k + 1. We need to show that for all a e NC n , X n0jb (a) = Y Cjh (a). Let 
l<m<2n — lbea number such that m <-^—> m + 1 . 
Case 1. Assume m ^ 2k + 2. 



Subcase la. Assume m $ b (that is, m ¥= bj for all j). Then by (33) 
we have that ® ai ,...,a k+1 ,b k+2 ,...,b n ( a ) = for an y a = ( a i> • • • > a fc+i) e A+i, so 

n,bW = o. 

Subsubcase la(i). Assume m is odd, i.e., m = 2j — 1 for some + 2 < 
j ^ n. If 7r e £ n (7To,b) then, since we assumed m $ b, 6j = 2 j — 2 and 
therefore by definition n m = —1, that is, ir(m) = 7r(2j — 1) < m, so in 
particular (m, m + 1) is not a little arc of it and therefore ^ n (cr) = 0. It 
follows that X^^ia) = = Y c ^{a). 

Subsubcase la(ii). Assume m is even, i.e., m = 2j — 2 for some k + 2 < 
j ^ n. If 7T 6 £ n (7To,b) then, since we assumed m $ b, = 2 j — 1 and 
therefore by definition 7r(m + 1) = vr(2j — 1) > m + 1 and in particular 
(m,m + 1) is not a little arc of tt. Again it follows that ^ n (a) = and 
therefore X no>h (a) = = Y ch (a). 

Subcase lb. Assume that m e b, i.e., m = bj for some k + 2 < j < n. 



Then by (33), for any a = (ai, . . . , a k +i) e *4fc+i we have 



$(a 1 ,...,a fe+1 ,6 fc+2 ,...,fe n )(^) =3™ $ (aibm) (0"m) 
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where b m is obtained from b by deleting the m and shifting all the numbers 
to the right of m down by 2. It follows that 

Y c , h (cr) = 3 n - 1 Y c ^Ja m ). (48) 

On the other hand, we have 

X 7ro , b {a)= J] J] 

n ire£ n (no,h), m«— ►m+l 

ire£ n (iTQ,h), m < / > m+1 

= 0+ ^ 3"- 1 * tm (a m ). 

7re£ n (7ro,b), m < - > m+l 

Relabeling 7r m in the last sum as it transforms it into 

7rs£' n _i(7ro,b m ) 



so by induction we get using (48) that X 7r0) b(c) = y^ 0i b(cr). 
Case 2. Assume m < 2fc + 1. 

Subcase 2(a). Assume m <— m + 1 (when thinking of 7To as "living" 
on the interval [2, 2fc + 1]), and note that in this case actually m must be 
< 2k. Then for any 7r e £ n (7T"o, b) we have 

^(a) = 3"- 1 ^ m (a m ), 

and consequently, denoting b' = (6^+2 — 2, . . . , b n — 2), we see that 

Xn ,b{°') = 3 n b ,(<7 m ). 

To evaluate Y c ^(a), we use ( |33| to get that 

K,b(0") = £ C a $ (a,b)(0-) = J] C a $ (a,b)(^) 

ae.A fe+1 ae_4 fc+1 ,msa 

= £ 3 n - 1 x(p(a,m))c a $ ( a m ,b')(^) 

ae^4fc_|_i, mea 

= £ da'$(a',b')(^m), (49) 



a'eA k 



where p(a, m) denotes the number of times m appears in a, a m is as in (32 ), 
and in the last line the sum was reorganized as a linear combination with 
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some new coefficients d a /. (The fact that this sum ranges over Ak explains 
why we needed to strengthen the inductive hypothesis at the beginning of 
the proof.) 

Note also that in a similar way, replacing the vector b with the empty 
vector of length 0, we have for any /i e NC^+i such that m * — * m + 1, 

Y cM = Yi C a $ a(/i) = J] C a $ a (/i) 
aeAk+i aeA k +i,mea 

2 3n_1 x(p(a,m))c a $ am (/x m ) 

ae*4i- +1 , mea 



2 d a >$ a >(fi r 



a'eA k 



with exactly the same coefficients d a i (the main thing to notice is that ap- 
pending the b at the end in the subscript of $ has no effect on the way the 



other hand, the case n = k + 1 discussed at the beginning of the proof gives 



that y Cj 0(/i) = X 7rf)i 0(/i) = 3 n l X- ( — r rxiP-m)- Combining this with (49), we 



recurrence (33) proceeds, i.e., the b "doesn't interact" with the a). On the 
other ha: 
that Y C) g 
get that 

a'eA k 

for any v e NCfc, or in other words simply that 

a'eAk 

Invoking the inductive hypothesis again, we conclude that 

a'eAk 

In particular, 

3 ( " _1) X^ mjb/ (0- m ) = J] da'$(a',b')Om)- 

a'e^lfc 

But as we showed above, the LHS is equal to X 7ro b (cr) and the RHS is equal 
to y Cj b(o"), so this proves the claim. 

Subcase 2(b). Assume m «-/-+ m + 1. Then for a similar reason as in 
subcase 2(a) above, X WOt b(a) = 0. We need to show that also Y Cj b(a) = 0. 
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The computation from subcase 2(a) is still valid, so we can write as before 
y c ,h{o) = 2a'G^ fc da'$(a',b')(o"m), and also 

Y c,0^) = Yj rf a'$a'(/i m ) 
a'eA k 

for any fi e NC^+i for which m * — > m + 1. But we know (again from the 
case n = k + 1) that Y c g{\i) = ^^(/j) = 0, so we deduce that 



2 rf a'$a' = 0. 



a'e.4 fe 

We claim that this implies also that 

da'$(a',b') = 0, (50) 



which, if true, would also give in particular that 



Yc,b(o) = d a /$( a / ib ')(Om) = 0, 

a'eAk 

and finish the proof. 

It remains to prove (50). First, for convenience relabel a' as a, b' as b and 
n — 1 as n. The claim then becomes that if b = (b^+i, ■ ■ ■ , b n ) satisfies bj e 
{2j — 2, 2j — 1} for all j and (d a )ae.4 fe are coefficients such that ^ aeAk d a & a = 
then Xiae^ifc ^a$(a,t>) = 0- The proof is by induction on n\ the idea as before 
is that the a' and b' components "don't interact" when we perform recursive 
computations using (33), and the proof involves similar ideas to those used 
above. If n = k there is nothing to prove, so assume that n > k. Let 
v e NC n , and pick some g such that g «-^-» g + 1. Divide into two cases, 
when g ^ 2k or when g < 2k — 1. We will show that in each case we 
have 2as^u^a < ^ ) (a,b)( zy ) = 0, which will imply the claim for the usual reason 
(Theorem 2.3). 

In the first case, if g $ {b k+1 , ...,b n } then Z! ae ^„ 4$(a,b)(^) = by Q. 
Alternatively, if g e {6^ +1 , . . . , b n } then, again by (33), we have that 



_ d a $ (a , h) (v) = T- L ^ ^<f (a ,b)(^), (51) 

aeAk aeAk 

where we denote b = (bk+i — 2, . . . , b n — 2) (note that the x(p) factors in 
(33) are all equal to 1 because of the structure of b). The right-hand side 
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of (51) is equal to by the inductive hypothesis (applied with the vector b 
replacing b). 

Finally, in the case where g < 2k — 1, we have 



2 4$(a,b)M 



aeA k 



2 cU"- 1 ^^))^^ 

ae.4 fe , ge{ai,...,a k } 



Vn 



(52) 



(where the notation p(a, (?) was defined immediately after (49)), and this can 
be rewritten as a sum of the form 



Yj /a' $ (a',b)(^ 
a'e«4 fc _i 



(53) 



where (/ a ')a / e^i fc _ 1 ar e some coefficients. On the other hand, by repeating the 
same computation in the case n = k where b is replaced by the empty vector 
of length 0, we see that for any fi e NC^ such that g «-^-» g + 1, we have 



2 = 2 /a'$a'(£ 

ae_4 fe a'eAk-i 



9/1 



(54) 



where the key observation is that the coefficients are the same numbers f a ' 
as in (53). But we know (by the assumption of the claim we are trying to 
prove) that the left-hand side of (54) is equal to 0. Since this is true for any 
ji 6 NCfc for which g * — » g + 1, we conclude that X/ a 'e.4 fc : /a'^a' = 0. The 
inductive hypothesis now implies that X!a'e.4 fe x /a'^a'E) = 0- I 11 particular, 
the expression in (53), which is equal to the left-hand side of (52), is 0. □ 



2.5 Proof of Theorem 1.13 



We are now in a position to prove Theorem |1.13[ with an explicit description 
of the multivariate polynomial F 7TQ (wi, . . . ,Wk) whose existence is claimed in 
the theorem. Let 7r e NCfc. The polynomial F no simply encodes the 7r th 
row of the matrix C n ; more precisely, we define it as 



Fn (wi, ...,w k ) 



=(ai,...,a fc )sNC+ 



k 

n 



C^al \W 3 



2j-dj 



(55) 



Denote f2 ra (z) = Y\i 



-Zi)(l+Zj+ZiZj). Taking the 1-evaluation of 



both sides of (44) and combining the resulting equation with Theorems 2.11 
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(l,l+Oi,...,l+Ofc)6fc+2,-j6n)- 



and |2.22[ we get that 

ASM(n)P(nW e ^(7r )) = £ A 

7rE£ n (7r ) 

a=(ai,...,a fc )eNC+ b k+2 ,...,b n 

b jS {2j-2,2j-l} 

The left-hand side is also equal to ASM(n)P (^n < II* n M by the rotational 
symmetry of II* n ' ) . To evaluate the right-hand side, represent each summand 



0(i,i+o l ,...,i+ofci6fc+2,-,6r») as a constant term using (37), to get that the last 
expression can be written as 



CT Z 



a=(ai,...,a fe )eNC + 



CT 7 



CT 7 



CT Z 



n ,T± 1 " 1 

6 fc+2 ,..,6„ ^ j ^ +2 ^- 

Vj 6 J -e{2i-2,2j-l} 

1 n / 1 1 

a=( ttl ,... A )ENC+ i=2 J j=k+2 \*j 

k+1 -. n 

2 ^Awn^ n 



a=(ai,...,a fc )eNC+ 



1 + Zj 

j=2 j=fc+2 



n n -, / fe+l 



2j'-2-a,_i 



j=fc+2 



J'=l J \a=(oi,...,a fc )6NC+ i= 2 



CT, 



2j-2 ' ^7ro(^2) • • • 5 #fc+l) 



j=k+2 



which is clearly equal to ASM(n) times the right-hand side of (15). □ 



Appendix A. Summary of algorithms 

The research described in this paper has been significantly aided by computer- 
aided experimentation and numerical computations performed by the author. 
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We believe use of such methods is likely to continue to play a role in the dis- 
covery of new results extending and building on this work. In the hope of 
stimulating such further research, we summarize the theory presented above 
in the form of explicit algorithms for computing quantities of interest re- 
lated to connectivity patterns of loop percolation and the theory of wheel 
polynomials. 

Algorithm A Computation of the probability vector /j, n = (yU 7r ) 7rS NC„ 
Step 1. Compute the matrix S n = (s 7r)0 -) 7riCre Nc, l where 

s*,o = #{1 ^ j < 2n : ej(ir) = a) 
(with ej defined by the equation (H), interpreted modulo 2n). 
Step 2. Compute /-i n as the solution to the vector equation 

H n (S n - 2nl) = 
(where I is the identity matrix), normalized to be a probability vector. 

Algorithm B Computation of the polynomials Q no (assuming Conjec- 



ture 1.12) 



Step 1. Use Algorithm A above to compute the sub matching event prob- 
abilities p n = P (tt < nt n) ) for n = k, k + 1, . . . , k(k + 3)/2. 

Step 2. Compute the numbers q n = p n 0^=i(^ 2 ~~ j 2 ) fc+1 ~-?. 

Step 3. Use the Lagrange interpolation formula to compute the unique 
polynomial G wo (m) of degree < k{k: + l)/2 satisfying 

G no (n 2 )=q n (k ^ n ^ k(k + 3)/2). 
Step 4. Q no is given by Q Wo (n) = G^n 2 ). 

Algorithm C Computation of submatching event probabilities in the half- 



planar model (assuming Conjecture 1.12) 
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Use Algorithm B above to compute the polynomial Q no . Let q* o be the 
leading coefficient of Q no . The probability F(n < U*) is g* o /2 fc(fc+1) . 

Algorithm D Computation of the matrices C n , C„ 

Translating (35) to a notation appropriate for the computation of C n)l7 , 
we denote Pj,h{^) = #{j ^ £ < k : 7r £ = 1} — — j — 1). The matrix 
C„ = (C , 7r,o-)7r,o-eNC„ is now computed using the explicit formula 



j<k, j^-^k 



(with x(-) defined in (34)), and C n = C^ 1 is its inverse matrix. 



Algorithm E Computation of the polynomials F^ (see Theorem 1.13) 
Step 1. Compute using Algorithm D above. 



Step 2. Use (55) to compute F 7 



7T0 ' 



Appendix B. Proof of Theorem 1.8 



We start by writing a more explicit formula for the entries of the transition 
matrices T„ . To do so, encode each row of plaquettes as a vector a e {0, l} 2n , 
where a coordinate corresponds to a "type 0" plaquette, defined as the 
plaquette shown on the left in Fig. |6j and a 1 coordinate corresponds to a 
"type 1" plaquette, which is the one on the right-hand side of the same figure. 
The parameter p corresponds to the probability of a type 1 plaquette. For a 
given plaquette-row vector a e {0, l} 2n , denote by / a (-) an operator that takes 
a noncrossing matching tt e NC n and returns a new noncrossing matching 
7r' = /a(7r) which is the result of composing the diagram of tt with the row of 
plaquettes, in an analogous manner to the composition of matching diagrams 
with Temperley-Lieb operators shown in Fig. |8j 

With this notation, it is clear from the definition of the transition matrix 
Tn that its entries (Tn ) 7rj7r ' are given by 

(Ti p \y = 2 p^^il-p)^ 1 -^ (7r,7r' e NC n ). 

ae{0,l} 2n 
/a(7r)=vr' 
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(a) 



(b) 



Figure 13: (a) Composing the diagram of a matching 7r with a row of type 
plaquettes causes a simple rotation; (b) flipping the kth plaquette from type 
to type 1 and leaving all other plaquettes as type is equivalent to the 
application of a Temperley-Lieb operator e k prior to the rotation. 



Differentiate this equation with respect to p at p = 0. On the right-hand 
side the nonzero contributions will come from the vector a = = (0, . . . , 0) 
and from the vectors of the form a = a& = (0, . . . , 0, 1, 0, . . . , 0) with only 
one coordinate in some position k equal to 1. This gives 

j- {T^\ y = -2n5 /o(7r)y + k < 2n : = tt'}. 

a P \p=0 



It is now not difficult to check (see Fig. 13) that fo(n) = p(n), where p 
is the rotation operator on matchings (defined in Theorem 2.7[ ), and that 
the operator / afc bears a simple relation to the Temperley-Lieb operator e k , 
namely, we have f ak = p o e k = e k _ x o p (where ek-i is interpreted with the 
usual mod 2n convention). 

Summarizing, since the vector pL n is the stationary probability vector for 
the matrices Tn \ that is, it satisfies n n Tn^ = U n , the matrix Q n = 4- Tn^ 

a P \p— o 

satisfies /-i n Q n = 0, and by the above observations it is easy to see that it is 
related to the operator H n from ^ by H n = ~R~ 1 Q n , where R is a matrix 
version of the rotation operator p, whose entries are given by R n y = 5 p ^)y. 
It remains to observe that pL n is invariant also under the rotation action, i.e., 
we have n n R = n n1 to conclude that li n H n = 0, as claimed. □ 
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